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ABSTRACT. It is shown that the shapes of acoustic radiation-induced static strain and 
displacement pulses (rectified acoustic pulses) are defined locally by the energy density of the 
generating waveform. Dispersive properties are introduced analytically by assuming that the rectified 
pulses are functionally dependent on a phase factor that includes both dispersive and nonlinear terms. 
The dispersion causes an evolutionary change in the shape of the energy density profile that leads to 
the generation of solitons experimentally observed in fused silica. 
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INTRODUCTION 

The propagation of a sinusoidal acoustic wave through a nonlinear medium is 
known to generate not only harmonics of the original waveform but also a static strain 
(acoustic rectification or “dc shift”) associated with the propagating wave. Using an 
approach based on the virial theorem, Cantrell [1,2] showed that the static strain 
generated by an acoustic toneburst (gated sinusoidal wave) is dependent on the energy 
density of the waveform. He considered that the energy dependence is local and 
predicted that the derived static displacement pulse would propagate as a right-triangular 
pulse as a direct result of the locality. Experimental confirmation of the local energy 
dependence of static displacement pulses was reported by Yost and Cantrell [3] and 
Cantrell et al. [4] for monocrystals of silicon for wave propagation along the pure mode 
propagation directions. Yost and Cantrell [3] also reported results for fused silica - a 
dispersive material with a negative nonlinearity parameter. The pronounced spreading 
and variations in amplitude of the displacement pulse profile observed in fused silica, 
however, were not predicted by the model. We develop here a model of static pulses in 
dispersive materials that does predict the spreading and amplitude variation of the 
displacement pulse profile as the consequence of the generation of solitary waveforms 
(solitons) of the acoustic energy density. 



ACOUSTIC RECTIFICATION IN NONDISPERIVE MEDIA 


We consider the propagation of an elastic wave in a lossless, non-dispersive, 
semi-infinite medium of arbitrary crystalline symmetry. The nonlinear equation of 
motion along an arbitrary propagation direction may be transfonned into the form 4 

( d 2 u / dt 2 ) = c 2 [ 1 - /3 s (du / da)](d 2 u / da 2 ) , where t is time, “a” is the coordinate along 
the wave propagation direction, u is the particle displacement, c £ is the wave phase 
velocity where s is a mode index associated with given directions of wave propagation 
and polarization, and [3 £ is the modal nonlinearity parameter. The Earnshaw particle 
velocity dujdt solution to the nonlinear wave equation for each mode s is [5] 

du 8 / dt = A s sin(<z>r - k s a + (f>) , (f) = k s a(/3 s / 2 c E ){du s I dt) (1) 

where ^4^ is the amplitude of the waveform, co is the angular frequency, k E = co/c £ , and (f) 
is a phase modulation factor, containing the nonlinearity parameter, that leads to wave 
distortion. For periodic waveforms of constant amplitude the acoustic strain rectification, 
or radiation-induced static strain <dujda>, is conveniently obtained by time-averaging the 
displacement gradient dujda (the angular brackets denote time-averaging). The strain 
rectification (static strain) is related to the time-averaged particle velocity <dujdt>, hence 
time-averaged energy density <E>, as [1,2] 

(dujda) = (/y 4 c 2 )((du £ /dt) 2 ) = (J./4jU,)(E) (2) 


where ju s is the modal second-order Huang coefficient. The last equality results from the 
relations E = pAduJdt) , where /y is the mass density of the propagation medium, and c £ 
= juJpo- 

It is clear from Eq.(2) that the static strain is defined by the time-averaged energy 
density <E> of the propagating wave. The constants /f. and /y serve as scaling factors for 
the static strain. Since the energy density for a toneburst in a non-dispersive, lossless 
medium propagates as a stationary waveform, the static strain must be local and therefore 
exists only in that region of space for which <E> is non-zero. Thus, from Eq.(l) the 

time-averaged energy density <E> = (l/2)po A z and the shape of the static strain pulse 
generated by a toneburst must be rectangular with height proportional to Ap. 
Accordingly, the static displacement is obtained as the spatial integral of the rectangular 
static strain pulse only over the non-zero values of <E>, or energy per unit mass A E , 
commensurate with the stationary wave nature of the propagating energy density (or 
energy per unit mass). The locality both of the static strain and displacement pulses 
means that, in contrast to harmonic generation, the static displacement does not increase 
linearly with total distance of wave propagation, but rather increases linearly only over 
the spatial extent of the non-zero energy density or energy per unit mass. Hence, for a 
toneburst of length L the static displacement pulse propagates in the shape of a right- 
triangular pulse with the peak value <u £ > obtained from Eqs.(l) and (2) as 
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FIGURE 1 . Rectified acoustic pulses generated by a 2ps toneburst: (a) Static strain pulse in non- 
dispersive media, (b) Static displacement pulse in non-dispersive media, (c) Static displacement pulse 
profile (soliton plus oscillatory tail) in dispersive media. 


(u £ ) = (\IS)P £ k 2 £ A 2 £ L . (3) 

Fig. la shows a typical rectangular static strain pulse generated by a 2ps toneburst 
as given by Eq.(2). Fig. lb shows a typical right-triangular static displacement pulse 
generated by a 2ps toneburst as given by Eq.(3). The perspective in Fig.l is the signal as 
seen by the receiving transducer. Yost and Cantrell [3] obtained experimental 
confirmation of the right-triangular shape of the static displacement pulse in a silicon 
monocrystal by substituting a right-triangular pulse from a function generator into the 
detection electronics of the experimental set-up and comparing the output signal with that 
generated from the acoustic wave. The overlap of the two output signals indicated that 
the static displacement was indeed right-triangular in shape. Calculations of the 
nonlinearity parameters from measurements of the slopes of the static displacement 
pulses along the pure mode propagation directions in single crystal Si gave /? E values that 
are in agreement with nonlinearity parameter measurements obtained independently from 
acoustic harmonic generation measurements [4]. Such agreement provides strong 
evidence of the right triangular shape of the static displacement pulses. 

ACOUSTIC RECTIFICATION IN DISPERSIVE MEDIA 

Further evidence that the static strain and displacement pulses are defined locally 
by the energy density of the propagating wave is obtained from a consideration of the 
static strain and displacement pulses in vitreous silica - a dispersive material with a 
negative nonlinearity parameter. The dispersion can be accounted by generalizing the 
phase factor (j) in Eq.(l) to include dispersive as well as nonlinear components. We 


assume that rectified pulses are generated by the toneburst in accordance with Eq.(2), 
where <E> = ( l/2)p 0 d, : 2 , to yield (ou / da) = {[) / 8c 2 )A }) , where we have dropped the 

modal subscripts and have replaced the energy per unit mass A by A n ■ For dispersive 
media the energy per unit mass A D 2 is no longer assumed to be constant but rather is 
assumed to evolve with propagation distance as A D = Av \koa- coot +(fiaj)}. At a = 0, we 
assume A r 2 = A 2 , the initial energy per unit mass of the toneburst launched into the 
material. The parameters co d and ko are associated with the toneburst and represent values 
of the parameters in the linear approximation. 

We assume that the nonlinear dispersion relation depends on A D and the wave 
number k as 


c o — co(k, A £) ) — (Oq + cq (k — k ^ ) + (1 / 2)cq ( k — k q ) 2 + Y Ajy + • • • , (4) 

where a Taylor series expansion with respect to ko and A if = 0 is used to obtain the last 
equality in Eq.(4). The Co, c'o , and T are expansion coefficients. We consider co and k to 
be dependent on time and propagation distance, and impose the geometrical optics 
postulate that co(a,t) = - ( dO/dt ) = ojq- (d<f>/dt) and k(a,t) = ( 86/da ) = k 0 + (d(f>lda), where 
6{a,t) = koti - ('X)t + (fAa,t). Substituting the expressions for oia,t) and k(a,t) into Eq.(4) 
and adding a spatial second derivative term to account for dispersion [6], we obtain 

(^/5t) + c 0 (^/5a) + (l/2)co(a^/aa) 2 +TA^ + (cq / 2A D )(d 2 A D / da 2 ) = 0 . (5) 

Transforming to a set of coordinates given by % = a-ctf and r = c(f , we re-write Eq.(5) 
in the form 

(d^/dr) + (l/2)co(5^/d£) 2 +(r/coM|) +(l/2A D )(d 2 A D /dg 2 ) = 0 . (6) 

Since we are interested in obtaining the propagation properties of the energy per 
2 2 
unit mass A D = E/p 0 , it is necessary to invoke an second equaton involving both A D ~ 

and(d^/ eg) . We thus invoke the conservation of energy expression 

dA 2 D ldt + d(c g A 2 D )lda = 0 (7) 

where from Eq.(4) c g = dco/dk = Co + ( d(f>lda ) is the energy propagation velocity (group 
velocity) in the linear approximation. Again transforming to the set of coordinates given 
by % = a-ctf and r = cfat , we re-write Eq.(7) in the fonn 

aA%/dr + (d/d€)[A%(d0/d$)\ = 0 . (8) 

Assuming that both A if and (dtp/ 8g) have stationary solutions to Eqs.(6) and (8) of the 
form / = f(%-V r) = f(z) where z = (E, - Vr) , we obtain from Eq. (8) 



(9) 


(d!dz)[VA 2 D - (d(j)ld£)A 2 D \ = 0 . 

2 2 

We integrate Eq.(9) and set the integration constant to VA~ where A is the energy per 
unit mass of the toneburst at a = 0. Solving the resulting expression for ( d(j)l 6c ) , we 
obtain 

£■ 

d<i>/dz = v( i + — ). (10) 

4 
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Following Karpman [7], we now assume that <j) = -// 0 z + (f)\(q - Vi) , substitute 
the expression into Eq.(6), and use Eq.(10) to obtain 

- nl - \ V 2 (1 - 4 V ) + (r / c' 0 )A 2 d + (1 / 2 A D )(d 2 A D /d{ 2 ) = 0 . (11) 

The stationary solution to Eq.(l 1) has the form of a solitary wave (soliton) given by [7] 

4 = 4 + 4 sec h 2 444 (4 - Vr) (12) 

where 

V = (T /c' 0 ) 1I2 (A$ + A 2 ) 112 (13) 

2 

and A s is the amplitude of the soliton. 

Strictly, Eq.(12) is the general stationary solution for a perturbation imposed on a 
continuous carrier wave of amplitude Aq and is valid for an arbitrary amplitude of 
perturbation [7]. If we assume that the perturbation (the toneburst) is large compared to 
the carrier amplitude, then in the limit where the carrier amplitude becomes arbitrarily 
small we obtain 


4 =A 2 sech 2 JT7^A s (t-VT) (14) 

1/2 

where V = (T/c'q) A s . The number of solitons generated, each having different 

amplitudes^ , depends on the length and amplitude of the initial perturbation. For a 
toneburst of amplitude A the initial rectification is a rectangular static pulse of 

amplitude A . The complete solution also includes an oscillatory tail wave packet similar 
in form to that generated by the Airy function. The static strain <du/da> is obtained 
directly from Eqs.(2) and (14) as <du/da> = <Gu/dg> = (fM/i)<E> where <E> = 
(l/2)poV. 

Fig.lc shows the static displacement pulse generated by a 2ps toneburst obtained 
by integrating Eq.(14) and the oscillatory tail with respect to £ The soliton is detected 



first, since it has the greater velocity, and is followed by the slower oscillatory tail. It is 
important to note that the rectified acoustic pulses are generated by the toneburst at a = 0. 
The dispersive and nonlinear properties of the rectified pulses are governed by the values 
of the expansion coefficients in the dispersion relation, Eq.(4), evaluated at the values of 
k 0 and (Oo associated with the toneburst. Thus, it should not be expected that launching a 
pulse in the material directly from a pulse generator would necessarily generate solitons, 
since the values of k (} and ox, used to evaluate the coefficients in Eq.(4) from such a pulse 
would not be the same as that of the rectified pulses generated by the toneburst. 

MEASUREMENT OF THE ACOUSTIC NONLINEARITY PARAMETER IN 
DISPERSIVE MEDIA 

Fig.lc shows the calculated static displacement profile, as observed by the 
receiving wave detector, resulting from a single soliton with an oscillatory tail generated 
by a 2ps tone-burst in a dispersive medium with [5 < 0. Note that since [S < 0. the 
displacement profile is inverted relative to that of Figs, la and lb where P > 0. The pulse 
profile of Fig.lc is also considerably longer than the toneburst duration and shows large 
variations in displacement amplitude. Both the spreading and the variations in amplitude 
of the static displacement pulse have been observed experimentally [3]. This observation 
provides further evidence that the static strain and displacement pulses are defined locally 
by the energy density of the waveform. 

We consider the implications of the static pulse profile to the measurement of the 
acoustic nonlinearity parameter in dispersive media. The initial portion of the static 
displacement pulse in Fig.lc results from the soliton contribution to the profile. The 
slope of the pulse profile following the first peak is given by <du/da> = (fi!4ju)<E>. At 
the half-peak height of that slope, where E, = Vx , the average energy density is evaluated 
from Eq.(9) to be <E> = (\/2)/x,A s . Generally, A s is not equal to the square of the 

amplitude of the tone-burst. However, for a single soliton generated by a toneburst of 

2 

appropriate length and amplitude it can be approximately so. If we assume that A s is 
very nearly equal to the square of the tone-burst amplitude, then measurement of the 
slope at f = Vi allows a direct, although approximate, calculation of fi. Yost and 

Cantrell [3] perfonned such a calculation of [i for Suprasil W1 vitreous silica, fortuitously 
at % = Vt, without benefit of the present analysis. Their calculated value of -12.7 for 

fi is in quite good agreement with the value -11.6 obtained independently from acoustic 
hannonic generation measurements for the material. The general agreement once again 
affirms the role of the local energy density in the generation of static pulses. 

CONCLUSION 

The analytical model previously developed for acoustic rectification in non- 
dispersive media [ 1 ,2] has been extended to the case of dispersive media. The extension 
is obtained by assuming that the energy per unit mass A d 2 of the rectified waves is no 
longer constant, as in the case of non-dispersive media, but now depends on phase factors 
similar to that obtained for the Eamshaw particle velocity solution to the nonlinear wave 
equation. We assume that at a = 0 (the spatial position at which the toneburst is launched 



2 2 * 

into the material) A [f = A (the energy per unit mass of the toneburst at a = 0). Assuming 
that the dispersion relation co (k, A d ) is a nonlinear function both of the wave number k 
and the rectified wave amplitude A D , we expand co (k, A n ) in a power series and use the 
geometrical optics postulate to obtain a Boussinesq type of equation. We add a 
dispersive term to the equation and solve for the energy per unit mass waveform after 
invoking the energy conservation equation. 

The solution predicts the generation of an energy per unit mass solitary waveform 
(soliton). According to the analytical model the static strain pulse profile is defined 
locally by the energy density of the toneburst and thus has the shape of the energy 
density, hence energy per unit mass, profile. The static displacement pulse profile is 
obtained by integrating over the spatial extent of the static strain profile to yield a profile 
that has been observed experimentally. The agreement between the pulse profile 
observed experimentally and the profile predicted by the model lends validity to the 
assertion that the acoustic radiation-induced static pulses are dependent locally on the 
energy density of the propagating waveform. 
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